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We extend a recently introduced semiclassical approach to calculating the influence of local lattice 
fluctuations on electronic properties of metals and metallic molecular crystals. The effective action 
of electrons in degenerate orbital states coupling to Jahn- Teller distortions is derived, employing 
dynamical mean- field theory and adiabatic expansions. We improve on previous numerical treat- 
ments of the semiclassical action and present for the simplifying Holstein model results for the finite 
temperature optical conductivity at electron-phonon coupling strengths from weak to strong. Sig- 
nificant transfer of spectral weight from high to low frequencies is obtained on isotope substitution 
in the Fermi-liquid to polaron crossover regime. 
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I. INTRODUCTION 

Despite many years of research, techniques to cal- 
culate the influence of strong vibronic interactions on 
electronic properties of metals and metallic molecu- 
lar crystals are urgently needed. Strong couplings of 
conduction electrons to local lattice fluctuations have 
been observed in certain phonon mediated superconduc- 
tors with high transition temperatures and manifest it- 
selves among other effects in unexpected mid-infrared 
transitions in the optical conductivity in these materi- 
als, like in Bai-xKxBiOs 1 and in the alkali-doped Ceo 
fullerides. 2 ' 3 Standard Migdal Eliashberg (ME) theory 4 
is not sufficient to explain these phenomena although typ- 
ical phonon frequencies of physical interest may remain 
small compared to typical electron energies. The rea- 
son for the failure is that ME theory assumes that the 
underlying electronic ground state can be described by 
Fermi liquid theory. However, if the electron-lattice cou- 
pling strength A = A/t exceeds a critical value (of or- 
der unity) the conduction electrons are believed to form 
"small polarons", 5 so that the electronic groundstate is 
fundamentally reconstructed. 

Signatures of polaronic behavior have been most accu- 
rately studied in the case of one or few electrons em- 
ploying the Lang-Firsov transformation. 6 However, its 
applicability to the physically relevant case of metallic 
densities is not established. A suitable method to study 
the interaction of conduction electrons with local lattice 
fluctuations may be the dynamical mean-field theory 7 
and in recent years an increasing number of studies have 
been published. In spite of the substantial simplifications 
which occur when the conduction electrons are described 
by local fluctuating fields the resulting equations are still 
very complicated. An exact solution is only possible in 
the case of a single polaron. 8 Numerically "exact" solu- 
tions are avaiable using quantum Monte Carlo methods 9 
in the antiadiabatic limit f2 ~ t of unclear physical rele- 



vance and using numerical renormalization techniques 10 
at zero temperature. Especially the latter is a promising 
approach at least for simplifying models like the Holstein 
model and may be extended in future to finite temper- 
atures. Semianalytical techniques have the advantage of 
being instructive and of being extensible to more realistc 
model systems of conduction electrons coupling to local 
lattice distortions including maybe even electron-electron 
interactions. These approaches 11-14 rely on a detailed 
analysis of the effective action based on the dynamical 
mean-field theory. 

Here, we extend a "semiclassical" approach 14 recently 
introduced by one of the authors to general models of 
electrons in degenerate orbital states coupling to Jahn- 
Teller distortions. The classical (static) phonon modes 
are treated exactly and the quantum (dynamical) modes 
are expanded to second order. The adiabatic expansion 
of the classical action allows for treating polaronic ef- 
fects in the physical relevant limit of small adiabatic 
parameter, from weak to strong electron-phonon cou- 
plings. However, the method is restricted to tempera- 
tures which exceed or are comparable to an energy scale 
set by a renormalized phonon frequency. Numerically, 
the method is implemented for the simplifying Holstein 
model of (spinless) electrons at half filling. We resolve 
previous numerical restrictions 1 connected with the an- 
alytic continuation of the local electronic Green function 
from the imaginary to the real frequency axis which al- 
lows for a more complete understanding of dynamical 
quantities, especially the optical conductivity. A treat- 
ment of coupling strengths beyond the classical polaronic 
instability becomes numerically tractable. Pioneering 
work in the calculation of the density of states and the 
optical conductivity at low temperatures has been done 
by Benedetti and Zeyher 12 and we provide a comparison 
to their results. The remainder of the paper is organized 
as follows. In Sec. II we apply the dynamical mean-field 
theory to models of metallic molecular crystals or met- 
als with strong Jahn- Teller distortions. The polaronic 
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instability for classical vibrations is discussed in Sec. Ill 
and known results in the opposite quantum limit at zero 
temperature are reviewed in Sec. IV. Section V gives 
the adiabatic expansion of the effective action and dis- 
cusses the numerical implementation to treat the Hol- 
stein model at half-filling, as well as the limitations of 
the approach. Numerical results for the density of states 
and the optical conductivity are presented and discussed 
extensively. We conclude in Sec. VI. 



II. MODELLING 

We study metallic molecular crystals whose Hamilto- 
nians are described by 



N 



Molecule 



E 



•f c 



C i,a C i.a 



(1) 



Here, the operator c\ a creates an electron in the molecu- 
lar orbital (MO) a on the i-th molecule. For simplicity we 
suppress the spin index and consider spinless electrons. 
However, our treatment can be easily extended to a para- 
magnetic state by including the appropriate spin degen- 
eracy. We will assume that the molecular orbital state 
are single particle like, e.g. provided by a Hartree-Fock 
treatment of the molecule. In particular, the hopping of 
electrons on and off a molecule shall not alter the local 
electronic level scheme. In this respect the assumption 
of spinless electrons helps to avoid unphysical large fluc- 
tuations in the local occupation number of the partially 
filled MOs giving rise to the conduction bands. 

In general, the high local symmetry of a molecule will 
lead to an orbital degeneracy of the interesting, partially 
filled MOs. Therefore -^Molecule describes the coupling 
between a /-fold degenerate electronic term and N Jahn- 
Teller active modes 



N 



^Molecule = \ J2 ( K Qi + M ^ 

( Wn(Qi,...,Q JV ) 



(2) 
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Here, the Qi denote the normal coordinates of the Jahn- 
Teller distortions. According to the Jahn- Teller theorem 
the leading term of the coupling constants Wij is linear 
in Q and different matrix elements are symmetry related. 
Employing group theory we may characterize different 
cases according to different irreducible representations of 
all possible point groups. In particular E®e denotes the 
coupling of a double degenerated electronic term with 
two Jahn- Teller modes. 

We do not attempt to solve (1) directly which is a very 
difficult problem. Instead we are analysing an effective 
action of the form 



^eff - ^ E E Qi("k)( K + Mwl)Qi(w_ k ) -nn (3) 
11 i=i k 

-Tr„ Tr orb In {c(i/„) - W\Q{v n - . 

This is the action of a single molecule embedded in an 
effective environment when the fermionic fields (which 
occur in quadratic form) are integrated out. The lat- 
ter give rise to the mean- fields c{v n ) which describe the 
conduction electrons and depend on odd Matsubara fre- 
quencies v n — 2TrT(n + 1/2). The Qi(uJk) arc bosonic 
fields (describing the Jahn- Teller modes) depending on 
even Matsubara frequencies uik = 2irTk. Note that the 
J/F-term in the logarithm is nondiagonal in both the Mat- 
subara and the orbital index. The action S e ff arises in a 
dynamical mean- field (DMFT) 7 description of the orig- 
inal lattice model Eq. (1) if we assume a structureless 
and fast fluctuating environment of a single molecule, 
i.e. if the local Green function can be averaged over the 
/ = n oib orbital degrees of freedom 
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^orb 



Tr orb G 



=loc 



(4) 



The orbital averaging was first introduced in Ref. 11. The 
DMFT treatment provides a self-consistency expression 
for the dynamic mean-fields c. The partition function 
may be written as a functional integral over the bosonic 
fields 



V[Q]cxp(-S cB ). 



from which the local Green function Q follows 



G{Vn) = 



S InZ 

Sc(u n ) ' 



(5) 



(6) 



The molecules are arranged on a three dimensional lat- 
tice. However, after averaging over the orbital degrees of 
freedom in (4) it may be sufficient to work with an ef- 
fective density of states of simple form, e.g. sem icircular 
po(e) = (1/AO Efe<^- e fe) - l/(27rf 2 )V4t 2 - e 2 . In this 
case the self-consistency equation (4) can be written in 
the form 



iuj n + fl 



t z 



71 orD 



■ Qn ({C n }) • 



(7) 



Eqs. (3), (5), (6) and (7) form a complete set of equa- 
tions which, in principle, can be solved for the orbital 
averaged local Green function on the imaginary axis of 
a molecular crystal with Jahn- Teller active local distor- 
tions. However, one cannot evaluate the path integrals 
exactly in Eq. (5) and further approximations are neces- 
sary. We adopt the adiabatic expansion of the effective 
action Eq. (3) which includes the leading term in the 
non diagonal part of the electron-phonon interaction as 
discussed in Sec. V. 
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The present modelling is not restricted to the de- 
scription of metallic molecular crystals as pioneering 
works on the "colossal magnetoresistance" manganites 
demonstrate. 11 However, when applying the present for- 
malism to others than molecular crystals one have to keep 
in mind that local distortions of the crystal unit cell do 
influence the neighboring cells. In particular, in transi- 
tion metal oxides E g symmetry distortions of the oxy- 
gen octahedra move oxygen atoms of several unit cells. 
Therefore, in this case the presented modelling is only 
applicable to phenomena which involve incoherent fea- 
tures of the electronic Green function, i.e. which involve 
high temperatures or high frequencies. 

III. CLASSICAL ACTION: ANALYTICAL 
INSIGHTS 

A good starting point to gain some insight in the 
physics described by the Eqs. (3) - (7) is to analyse the 
effective action Eq. (3) in the high temperature limit 11 
when only thermal fluctuations are present. S e s reads 

K N 

S c i ass (Q) = — ^Q 2 - M n ( 8 ) 
1=1 

-Tr„Tr orb In [c{v n ) - W\Q]} 

where the Qi = Qi{tOk = 0) are the classical components 
of the bosonic fields. The function 

P(Q) = exp[-S* class (g)] (9) 

is the (unrenormalized) probability that the classic Jahn- 
Teller distortion takes the value Q. 

We specialize now to different cases. The case of a non 
degenerate electronic term coupling to a totally symmet- 
ric distortion (Holstein model) is the simpliest situation 
with the interaction matrix being a scalar W_ = g Q. The 
strength of the electron-phonon coupling is measured by 
A = g 2 / K. (In the following we will frequently denote the 
dimensionless coupling parameter by A — A/t). It is easy 
to show 11 that the effective action S c i ass (Q) changes from 
a single to a double well behavior with increasing A as 
shown in Fig.l. In Fig.l the term — Q/2 has been added 
to the effective action to symmetrize the distortions in 
the case of hole and electron occupation on a given site. 
A double well occurs when the <3 2 -cocfficicnt in the ef- 
fective action vanishes. This is the case for A c = — l/r 2 
with 

r 2 = T^ c - 2 K). (io) 

n 

Here the fully self-consistent determined dynamical fields 
have to be used. If we insert the free form, the critical 
electron-phonon coupling reads 15 : 

A c = ^ [1 - (M/2i) 2 ]" 3/2 • (11) 



High temperatures (J/Q. > 1) 




Q 



FIG. 1. MAIN PANEL: The (classical) effective action S cH as 
function of the distortion Q for the Holstein model at half-filling. 
Shown is the behavior for different coupling strengths A at a high 
temperature of T/t = 0.2. The forming of a double well at a 
finite value of A is clearly seen. INSET: Adiabatic potential energy 
surface Eapes(Q) = Q 2 /2A+Q(n — 1/2) at A = 2.56 for a molecule 
with no and one electron in the relevant molecule orbital. The 
distortion equals ±A. 

At couplings exceeding A c the motion of an electron is 
accompanied by a lattice distortion, i.e. small polarons 
are formed. It is important to note that the behavior of 
an isolated molecule is different. With a linear electron- 
phonon coupling a finite Jahn- Teller distortion occurs for 
every finite coupling constant A as illustrated in the inset 
of Fig.l. In a crystal the quantum fluctuations of the 
occupation numbers of the MOs due to the hopping of 
the electrons smear out this distortion. A finite minimal 
coupling strength is required for the lattice relaxing to 
the actual occupation of a given site. 

The occurence of finite Jahn- Teller distortions at a fi- 
nite critical value of A is a generic feature of the clas- 
sical action S c iass- However, the critical value depends 
upon details of the Jahn- Teller coupling, in particular 
upon the degeneracy of the problem. This can be illus- 
trated by comparing the T\ ® /i-problem describing the 
coupling between electrons and vibrations of the charged 
Ceo molecules to the E ® e-problem relevant for transi- 
tion metal oxides with perovskite structure. Electrons 
added to an isolated fullcrcnc molecule belong to the T\ u 
lowest unoccupied MO which can couple linearly to vi- 
brations of symmetry h g . 16 The Jahn- Teller interaction 
matrix reads 

/Qi-VSQa -V3Q 3 -V3Q 2 \ 
K=l\ -V3Q 3 Qi + VSQi -V3Q 5 • (12) 
V -V3Q 2 -V3Q 5 -2Qi / 

When introducing a new parametrization and expressing 
Qi,....Qs via one amplitude Q and four angles 16 the 
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Ql-type 



Q2-type 



Q5-type 








FIG. 2. The five /i-type Jahn- Teller distortions measured rela- 
tive to a sphere. 

classical effective action can be written in compact form 



Q2 

Soft = — — — — Tr„ In 
2AT 



{vn) ~ j Q c{v n 



Q 3 

— — cos 3a 
4 



(13) 



Expanding S e g for small Q reveals that the Q 2 -coefficient 
of the classical action vanishes for 



2 1 

3 fl 



(14) 



For coupling strengths A > A c finite Jahn- Teller dis- 
tortions occur as displayed in Fig. 2. For the E eg) e- 
problem one obtains 11 A c = — 1/2F2. Wee see that for 
comparable values of the "electron bubbles" T2 the oc- 
curence of finite Jahn- Teller distortions require stronger 
electron-phonon couplings in the higher degenerate case. 
In this sense the weak-coupling regime for the T\ ® h- 
problem extends to larger coupling strengths. 

It is important to note that the occurence of a po- 
laronic instability indicates a break down of Migdal- 
Eliashberg theory because the ground state is fundamen- 
tally reconstructed. Physically speaking a quasiparticle 
picture does not hold in the transition regime because 
neither a description as electrons nor as polarons is valid. 
This give rise to interesting isotope effects on electronic 
properties as discussed in Ref. 14 and in Sec. V. 



IV. INFLUENCE OF QUANTUM 
FLUCTUATIONS 



Quantum fluctuations of the lattice will change the 
value of the coupling constant A c needed to form small 
polarons. Even for stronger electron phonon couplings 
when a classical lattice would adjust to the occupancy 
of a given site, possible tunneling may hinder the forma- 
tion of polarons or even prevent it. At sufficiently low 
temperatures only zero-point quantum vibrations of the 
lattice are present and the classical picture is not valid. 
Detailed studies exist 10,12 ' 15 for the evolution of the elec- 
tronic effective mass in the Holstein model as function of 



7=0.2 (NRG) 
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FIG. 3. Inverse ratio of the effective and bare electronic mass 
as function of A. The classical prediction is compared to the results 
obtained from an NRG study by Meyer, Hewson and Bulla in Ref. 
10 at finite phonon frequency. In the shaded parameter regime 
strong isotope effects have to be expected. 

the coupling strength A in the adiabatic regime. All re- 
sults indicate that the electronic effective mass diverges 
(or at least increases very rapidly) at a finite value of 
the electron phonon coupling even in the presence of 
zero-point fluctuations. In Fig. 3 we replot the re- 
sults of a numerical rcnormalization (NRG) study pub- 
lished in Ref. 10. The inverse ratio of the effective and 
bare electronic mass as function of A (= 2g 2 /ujo in the 
notations of Ref. 10) is shown. For the adiabatic pa- 
rameter 7 — fl/t — 0.2 (t = 0.25) a critical value of 
A c (7) w 1.57 has been reported. The classical approach 
predicts m* /m = 1+A p(n) with A = A/(l — A/A c ), i.e. a 
diverging electron mass at the classical value A c which 
reads for a semicircular density of states at half filling 
A c = 37rf/[(2s + 1)4]. The classical result is also shown 
in Fig. 3. Note that in the NRG investigations the case 
of twofold spin degenracy is considered, i.e. the classical 
polaronic instability occurs at A c « 1.18. We conjec- 
ture that with increasing values of the phonon frequency 
the "polaronic instability" is shifted to larger values of 
the electron phonon coupling and in the shaded param- 
eter regime of Fig. 3 strong isotope effects on electronic 
properties have to be expected. 

If one includes the leading frequency correction to the 
finite loop (ME) diagram, the effective mass reads 15 



m 
m 



1 



l + Ap(/i)--A7 



(15) 



with 



7 being the renormalized adiabatic parameter 
7 = 7 y/l — A/Ac. Therefore, the electronic effective 
mass is basically unchanged for phonon frequencies in 
the adiabatic regime. Systematic adiabatic (small-7) 
expansions 15 at zero temperature give all contributions 
of order 7 and correct Eq. (15). However, these contri- 
butions come with additional factors A" (n > 2). Thus 
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they are not small close to the polaronic instability and 
an adiabatic expansion is not possible. This is one rea- 
son why the semiclassical approach discussed here is re- 
stricted to high to intermediate temperatures which help 
to re-establish smaller lattice distortions. 

The results of Benedetti and Zeyher 12 (although at fi- 
nite temperatures) for 7 = 0.01 are close to the classical 
result and fall as they should in the shaded region of Fig. 
3. One of the authors has studied the consequences of 
the isotope shift of A c when both thermal and quantum 
fluctuations are competing by means of the semiclassical 
approach revisited here. As one of the most dramatic 
manifestations an isotope driven insulator-to- metal tran- 
sition may be observed. 14 



V. ADIABATIC CORRECTIONS 

From the above considerations it becomes obvious that 
adiabatic corrections to the classical effective action have 
to be taken into account. For the Holstein model this has 
been done in Ref. 14. The idea is to expand the difference 
S'off-S'ciass to second order in (l—6 nm ) W-{v n -v m ), i.e. to 
leading order in the non diagonal part of the electron- 
phonon interaction. This yields 



Intermediate Temperatures (J/Sl ~ 112%) 



Soft — ^class + S2 



(16) 



with 



^2 = ^ E QM(K + Mu>l)$(u>- k ) + (17) 



k^O 



x [c(u n - k )-m$j\ W\Q{-u k )\ 

In the following we will present numerical results for 
the Holstein model at half-filling, i.e. W_ = gQ and fj, = 0. 
We generalize the treatment given in Ref. 14 by using a 
different numerical implementation which does not rely 
on approximatively continuing the local electronic Green 
function from the imaginary to the real frequency axis. 
Instead, the continuation is done "numerically" exact. 
This allows for a more extensive discussion of the opti- 
cal conductivity. We shortly review results from Ref. 14 
pertinent to the present work. For W_ — g Q all but one 
integration in Eq. (5) can be performed analytically and 
the effective action depends only on the classical coordi- 
nate Q 



5*cff(g)-5'clas S (0)+E ln 



fe>0 



1 



(18) 



Here, not the bare phonon frequency f2 2 = K/M enters 
but a renormalizcd one which depends on both the clas- 
sical distortion Q and the even Matsubara frequencies 
uo k = 2nT k and reads 



O 




FIG. 4. The effective action S e g as function of the distortion 
Q for the Holstein model at half-filling. The adiabatic parame- 
ter equals 7 = 0.2. Shown is the behavior for different coupling 
strengths A at a intermediate temperature of T/t = 0.036. The 
forming of a triple (instead of a double) well at a finite value of A 
is clearly seen. 



n 2 (uj k ,Q) = n 2 [i + n(w fc ,Q)] 



(19) 



The (jj k dependence arises from the fact that the phonon 
propagator is rcnormalized by electron bubbles, i.e. pro- 
cesses where one electron absorbs (or emits) one phonon 
and gains (or loose) the energy w k . Therefore in Eq. (19) 



n( Wfc ,Q) = AT^ 



c{v n ) - Q c(v n -u>k)-Q 



(20) 



If the renormalizcd phonon frequency were independent 
of the even Matsubara frequencies to k — 2nTk, the quan- 
tum term S2 were the action of a harmonic oscillator. 
This simplification was explored in Ref. 14. Here, how- 
ever, we will investigate the exact form of the second 
order quantum corrections. 

At weak couplings the mean-field functions c{vn) may 
be approximated by the free ones. However, at interme- 
diate to strong electron phonon couplings the electronic 
properties are influenced and the mean-field functions are 
renormalizcd as described by the self-consistency equa- 
tion (7). The local Green function is obtained by the 
logarithmic derivative of the partition function Z with re- 
spect to c(v n ), Eq. (6). It reads (denoting i J dr P(Q) 

by(» 



Q{Vn) 



1 +F{v n ,Q) ' 



c{y n ) - Q 



c{v n ) - Q , 



(21) 



Comparing to the classical result, the quantum correc- 
tions enter in two places. First, they enter in the modified 
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probabilty P(Q) = cxp(— S c g) that the classic distortion 
takes the value Q. Second, they give rise to an additional 
scattering term F(v n , Q) 



F(u n ,Q) = ATj2 



n 2 



k>0 



wl + ni(Q) 



(22) 



1 



+ 



1 



c{v n -k) - Q c{u n+k ) - Q 



which can be interpreted as generalized Migdal self en- 
ergy where the classical k = term is subtracted. The 
two modifications have opposing influences. First, as 
shown in Fig. 4, the adiabatic corrections lead to an in- 
creased propability for distortions around Q — com- 
pared to the classical one. In consequence, quantum 
fluctuations of the lattice (not present in the classical 
approach) may lead to an increased mobility of the con- 
duction electrons. Second, when an analytic continua- 
tion to the real frequency axis of F is performed [yield- 
ing F(lo,Q)], a finite imaginary part implies additional 
electron phonon scattering. Note, that for small frequen- 
cies lo below the renormalized one the imaginary part of 
F(lo, Q) may be reduced or vanishing. 



A. Numerical implementation 

Iteration of Eq. (7) and Eq. (21) will yield the local 
electronic Green function Q(v n ) on the imaginary axis. 
The occuring fermionic and bosonic Matsubara sums are 
calculated introducing high frequency cut offs Co and v, 
respectively. The high frequency contributions to the 
sums are approximated by integrals. For reference we 
summarize here the results where functions f(Co) {f{v)} 
are finite bosonic (fermionic) Matsubara sums up to uo k < 
Co (up to v. n < v): 



£L 2fi Co 
S 2 (Q) = S 2 (Q,ui) + — - — arctan - 

co , / n 2 

In 1 + 



2ttT 

U(uj k ,Q) = U(co k ,Q,D 



F{v n ,Q) = F(v n ,Q,Co) 



2nu; k V Q 2 + D 2 

An i 



7T Q 2 -n 2 



(23) 



(24) 



(25) 



LO IT 

Q (arctan - - - ) 

_ , Cd TT . 

il (arctan — — sgnt^ — J 



In principle, the truncation at Co and the adding of an 
asymptotic contribution in Eq. (25) is not valid if v n w Co, 
i.e. for large values of the fermionic Matsubara frequency. 
However, the quantum scattering term F(u n , Q) tends to 
zero in this case and small errors are not important. In all 



our calculations we used 128 bosonic (fermionic) Matsub- 
ara frequencies. To speed up the integration performed 
in Eq. (21) the electron electron bubble U(uok,Q) has 
been interpolated as function of Q for each value of uJk 
in each DMFT iteration step. 

In the present approach only one dimensional integrals 
over a real variable have to be calculated. This can be 
done with high precision yielding Q{v n ) up to 14 digits. 
Therefore, we can continue the electronic Green function 
to the real frequency axis employing a "simple" Pade 
approximation. 17 The quality of the results are even more 
improved when using the classical solution as a refer- 
ence state, i.e. doing the continuation for the difference 
G — ^ciass- This is a real strength of the present semi- 
classical approach allowing the calculation of dynamic 
quantities without further assumptions. As important 
example we have calculated here the finite temperature 
optical conductivity of the Holstein model as convolution 
of two full Green functions: 



cr(Q) = (T 



x J depo(e) p(e, uo) p(e,to + 0) 



(26) 



The external frequency fl shall not be confused with the 
phonon frequency. In Eq. (26) is po{e) of semicircular 
form and the spectral function defined by 



1 



1 



LO + H - S(w) - 



(27) 



The self-energy S(w) can be obtained from the local elec- 
tronic Green function G(u>) on the real frequency axis. 
We also present results for the density of states 



p(io) = J dep (e)p(e 



LO) 



(28) 



To summarize, our implementation of the scmicalssical 
approach for the Holstein model differs in two way from 
the one discussed in Ref. 14. First, when calculating the 
effective action Eq. (17) we do not rely on a static ap- 
proximation, i.e. we do not replace Cl(u>k,Q) 0(0, Q). 
However, this has little influence on the results. Second 
and most important, G{v n ) is continued to the real axis 
not approximatively but numerically exact. In particu- 
lar, the investigations of dynamic quantities for electron- 
phonon couplings above the (classical) polaronic insta- 
bility is possible which is not achieved in Ref. 14. 



B. Limitations 

Before presenting results we shortly summarize the 
shortcommings of the presented semiclassical approach as 
discussed extensively in Ref. 14. Within the framework 
of dynamical mean-field theory the approach becomes ex- 
act at high temperatures T — > oo, irrespectively which 
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value the electron phonon coupling takes. However, at 
low temperatures the approach fails for two reasons: 

First, at intermediate to strong couplings the weak 
coupling (Fermi liquid) ground state is fundamentally 
reconstructed. In the present formalism this is re- 
flected by the occurence of imaginary phonon frequen- 
cies f2 2 (wfe,Q)- In this case 52 (<2) is not well defined 
anymore at low temperatures. Only larger temperatures 
re-stabilize the ground state. To explore this in more de- 
tail, let us suppose the renormalized phonon frequency 
\fl(u>k = 0, Q = 0)| turns imaginary at a critical coupling 
A c . For T — > and for couplings just above the critical 
value AA = A - A c one finds ft 2 (0,0) = -O 2 AA/A C . 14 
Because tt 2 (u)k,Q) is a strictly monotonic decreasing 
function with Matsubara frequency and displacement, 
the quantum correction 52 (Q) is only defined if 



(29) 



However, at high temperatures ft 2 (0,0) = ft 2 - A/(4T) 
and the renormalized phonon frequency is always real. 

Second, also at weak couplings the ground state cannot 
be reached. A purely classical approximation (fl — > 0) 
would lead to a T-linear scattering rate with coefficient 
of order unity. The semiclassical approximation overcor- 
rects for this behavior, cancelling the leading term leav- 
ing a correction of order 7T but of the wrong sign. Note 
that for T < fl this is 0(-/ 2 ). Note also that at T = 
the adiabatic expansion is possbile for weak couplings 
(see Sec. IV and Ref. 15). 

In conclusion the semicalssical approach is restricted to 
a range of high (77 £1 l/27r) to intermediate (T/Cl <~ 
l/27r) temperatures and we present only results in this 
range. 



C. Results 

In the following we present numerical results for the 
density of states (DOS) and the optical conductivity at 
finite temperatures for electron-phonon couplings from 
weak to strong. The model under consideration is the 
Holstein model with one (spinless) electron per each two 
sites (half filling), i.e. we discuss the simpliest case of non- 
degenerate electrons coupling linearly to non-degenerate 
lattice distortions. 

In Fig. 5 we show typical results at weak couplings, 
i.e. for values of A well below the (classical) polaronic 
instability of A c w 2.36. We recover the behavior known 
from calculations employing Migdal-Eliashbcrg theory. 18 
A tiny quasi-particle peak at zero frequency (and the in- 
set in the upper panel showing the DOS over the complete 
frequency range shall demonstrate that it is indeed a tiny 
structure) is observed in the DOS and no upper and lower 
sub-bands are formed. The quasi-particle peak originates 
from the fact that conductions electrons can only inelasti- 
cally scatter from quantum lattice fluctuations with a fi- 



Weak coupling 




FIG. 5. UPPER PANEL: DOS of the Holstein model 
(half-filling) at a weak electron-phonon coupling A = 1 and a tem- 
perature T/t = 0.04. Results for various values of the adiabatic 
parameter 7 are shown. For finite phonon frequencies a small 
quasi-particle peak is visible at u> = 0. The inset displays the 
DOS over the full frequency range for the largest value of the adi- 
abatic parameter, 7 = 0.1. On this scale the quasiparticle peak is 
hardly seen. LOWER PANEL: The frequency dependend optical 
conductivity for the same model and parameters as in the upper 
panel. It displays a Drude-like behavior. 
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nite frequency transfer. Therefore, the peak gets slightly 
broader with increasing adiabatic parameter 7. Note that 
in all our calculations 7 is a small parameter. Following 
the behavior of the DOS, the optical conductivity shows 
a typical metallic behavior with a sharp Drude peak. The 
optical conductivity drops monotonically with increasing 
applied frequency u. 

Fig. 6 shows the DOS and the optical conductivity at 
an intermediate electron-phonon coupling, just above the 
(classical) polaronic instability. In this regime isotope ef- 
fects (i.e. changing values of 7) may alter the behavior 
qualitatively. At low phonon frequencies sub-bands in 
the DOS are weakly developed. With increasing 7 these 
sub-bands disappear in favor of a quasi-particle peak at 
uj = 0. The quasi-particle peak is more pronounced as 
in the weak coupling case with slightly smaller width. 
With increasing adiabatic parameter also the behavior 
of the optical conductivity changes qualitatively from 
an insulating behavior with increased intra-band absorp- 
tion to a Drudc-likc metallic behavior with an incoher- 
ent contribution visible as "shoulder" . Although a Drudc 
peak develops the conductivity is one order of magnitude 
smaller than for small electron-phonon couplings. One 
may speak of a "bad metal" . 

Results for stronger couplings are shown in Fig. 7. 
Sub-bands are well developed in both the DOS and the 
optical conductivity. Only for a significant change in the 
phonon frequency a Drude-like metallic behavior is re- 
stored. In this sense isotope effects are less pronounced 
in the strong coupling case as in the intermediate one as 
they do not lead to immediate qualitative changes. The 
qualtitative behavior as given in the classical picture is 
maintained for all small values of the adiabatic param- 
eter 7. However, at larger phonon frequencies both a 
quasi-particle peak at oj = and sub-bands are visible 
simultanously, a feature which cannot be obtained within 
Migdal-Eliashberg theory. The intra-band absorption in 
the optical conductivity occurs around u) ~ A and shifts 
to smaller frequencies with increasing value of the adi- 
abatic parameter. The intra-band absorption is asym- 
metric about its peak. In particular, the low-frequency 
side of the absorption peak is elevated above the high- 
frequency side. 

The low temperature behavior of the DOS and the 
optical conductivity in the adiabatic regime has been 
carefully studied by Benedetti and Zeyher. 12 Their re- 
sults confirm nicely that our findings at intermediate 
temperatures are not spoiled by the low temperature 
break down of the semiclassical approach but that they 
reflect the physical behavior of the model. In partic- 
ular, the optical conductivity shows the behavior pre- 
sented here: (i) Drude like behavior at weak coupling, 
(ii) quasi-particle peak together with incoherent features 
at intermediate couplings and (iii) intra-band absorp- 
tion at strong coupling. Obviously, the present semi- 
classical approach is capable of describing the interesting 
crossover regime from high to low temperatures. Un- 
fortunately, isotope effects on the optical conductivity 
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FIG. 6. UPPER PANEL: DOS of the Holstcin model 
(half-filling) at an intermediate electron-phonon coupling A = 2.89 
and a temperature T/t = 0.04. Results for various values of the adi- 
abatic parameter 7 are shown. Increasing the adiabatic parameter 
has a considerable influence. Weakly developed sub-bands occur- 
ing at small phonon frequencies disappear in favour of a narrow 
quasi-particle peak at u> = with increasing 7. LOWER PANEL: 
The frequency dependend optical conductivity for the same model 
and parameters as in the upper panel. At small phonon frequencies 
an intra-band absorption is clearly visible. With increasing phonon 
frequency this feature merges with a developing Drude peak giving 
rise to a "shoulder" . 
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Strong coupling 
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FIG. 7. UPPER PANEL: DOS of the Holstein model 
(half-filling) at a strong electron-phonon coupling A = 3.24 
and a temperature T/t = 0.04. Results for various values 
of the adiabatic parameter 7 are shown. For all values of 7 
the DOS displays well developed sub-bands. A quasi-particle 
peak becomes visible for larger values of the phonon fre- 
quency. LOWER PANEL: The frequency dependend optical 
conductivity for the same model and parameters as in the 
upper panel. An asymmetric inter-band absorption is clearly 
seen situated roughly at ui ^ A, shifting towards smaller fre- 
quencies with increasing 7. 



were not discussed in Ref. 12 which are especially rele- 
vant when Migdal-Eliashbcrg theory breaks down (i.e. at 
intermediate electron-phonon coupling strengths). 

It should be noted that the qualitative features and 
changes in the optical conductivity when increasing the 
electron-phonon coupling from weak to strong are simi- 
lar to the ones predicted by a small-polaron absorption 
theory. 19 However, in this theory the hopping of the con- 
duction electrons is a small parameter contrary to the 
approach discussed here which is based on an adiabatic 
expansion. Therefore a comparison is only meaningfull in 
the limit of large electron-phonon couplings where small 
polarons are well defined quasi-particlcs. Here, the main 
absorption process is the excitation of a self-trapped elec- 
tron from its localized state at one site to a localized state 
at a neighboring one. Due to thermal and quantum fluc- 
tuations of the lattice the energies of the localized states 
are broadened and one obtains for the optical conductiv- 
ity 



o ~ — exp[-(A - ujf/U AE vih ) 



(30) 



where — 0/2 is the zero-point vibrational energy 
at low temperatures and i? v ib = T in the case of high 
enough temperatures when the vibrations can be treated 
classically. The classical (high temperature) result can be 
compared to our numerical datas in the strong coupling 
limit and reasonable agreement is found. 

Depending on the values of the phonon frequency f2 
and coupling strength A qualitative features of various 
experiments are present. To name two examples: (i) In 
Bai_ x K x Bi03 for high enough doping x (so that charge 
ordering is destroyed) the optical conductivity displays 
a Drude-like peak with a pronounced incoherent feature 
in the mid infrared; 1 (ii) Asymmetric intra-band absorp- 
tion can be observed in the alkali doped fulleridcs K3C60 
and RbaCgo- 2 As outlined in this paper the semiclassi- 
cal approach can be extended to more realistic models of 
electrons in degenerate orbital states coupling to Jahn- 
Teller distortions allowing for a more detailed analysis of 
experiments. Along the lines of Refs. 13 and 20 effects of 
local Coulomb repulsion may be included. 



VI. CONCLUSION 

In this paper we have extended a semiclassical 
approach 14 based on the dynamical mean-field theory 
recently introduced by one of the authors to treat the 
coupling of conduction electrons to local Jahn- Teller dis- 
tortions in metals and metallic molecular crystals. We 
have improved previous numerical treatments allowing 
a detailed discussion of dynamic quantities such as the 
finite temperature optical conductivity and presented re- 
sults for the simplifying Holstein model. In particular, 
the density of states and the optical conductivity could 
be studied not only for electron-phonon couplings below 
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the (classical) polaronic instability but also above. At 
intermediate to strong couplings the optical conductiv- 
ity displays features which are known from experiments 
on phonon mediated superconductors with high transi- 
tion temperatures. Depending on paramters, a Drudc- 
like peak with a pronounced incoherent "shoulder" or an 
asymmetric intra-band absorption is observed. Strong 
isotope effects occur at couplings close to the polaronic 
instability changing the qualitative behavior of the sys- 
tem. 

The presented results demonstrate that the analysis of 
simplified effective actions based on the dynamical mean- 
field theory and adiabatic expansions is a powerful tool 
to study the influence of local lattice fluctuations on elec- 
tronic properties of metals. Future work will reveal if the 
presented approach is capable also of a quantitative anal- 
ysis of electronic properties of metals and metallic molec- 
ular crystals strongly influenced by local Jahn- Teller dis- 
tortions. 
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